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Abstract. We demonstrate that two remote qubits can be entangled through an optically 
active intermediary even if the coupling strengths between mediator and qubits are different. 
This is true for a broad class of interactions. We consider two contrasting scenarios. First, we 
extend the analysis of a previously studied gate operation which relies on pulsed, dynamical 
control of the optical state and which may be performed quickly. We show that remote 
spins can be entangled in this case even when the intermediary coupling strengths are 
unequal. Second, we propose an alternative adiabatic control procedure, and find that the 
system requirements become even less restrictive in this case. The scheme could be tested 
immediately in a range of systems including molecules, quantum dots, or defects in crystals. 
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1. Introduction 

The control and physical representation of entanglement lie at the heart of quantum 
computing. This leads to several design considerations. First, it is essential to have a well 
defined localized information carrier that can be manipulated precisely. The basic information 
unit is usually a qubit - a two level quantum system. Second, it is necessary to design 
quantum gates to precisely control both the individual dynamics of each qubit, and their 
correlated motion. These two classes of operation are usually treated separately, and define 
single and two qubit gates. Together they form a universal set that can encode any quantum 
algorithm. Any practical quantum computer must, however, provide significantly more than 
just a universal gate set. It must be capable of preparing an initial quantum state, and of 
reading a final state. It should be able to control sufficiently many qubits for a time long 
enough to perform a useful calculation, without significant loss of entanglement during the 
evolution of the quantum state. And, as any quantum device will surely be controlled by a 
classical computer system, compatibility with current technologies is desirable. 

These requirements have led to a wealth of research on solid state implementations of 
quantum computers. In particular, an electron or nucleus with spin 1/2 constitutes a perfectly 
defined qubit. However, in solid state systems there is often uncertainty about the positions 
of and interactions between such spins, which of course makes control difficult. In this paper 
we shall demonstrate that it is nonetheless perfectly possible to devise very accurate quantum 
gates even when the parameters describing the qubits have a degree of randomness. 

Many papers have discussed how spins could be used as qubits in various materials 03 El- 
It is often the case that spins with the long coherence times that are so desirable for storing 
quantum information, are hard to control directly. However, they can be accessed indirectly 
through other states with shorter lifetimes, and manipulated more quickly. For example, 
different defects in semiconductors can be chosen such that have the correct properties play 
different roles: A defect with an electron spin in a quiet environment can be used to represent 
that quantum information, and interactions can be provided by optically active defects with 
shorter decoherence times 0. Alternatively, excitons in quantum dot systems can be excited 
to provide coupling between dot based spin qubits [@][51[6]|, that can sometimes be mediated 
by an optical cavity J71[8l. Further, in NMR quantum computing, the interaction between 
nuclear spins is provided by the electrons iH OH. A summary of such approaches can be 
found in Table 1. 

More recently, several key experiments have been performed that demonstrate many of 
the ingredients that are needed for the operation of these control schemes. For example, the 
spin of an NV- centre in diamond can be initialized, read and manipulated optically [11], 
and can be coupled to other nearby electron spins coherently, and this coherence can be 
manipulated optically [[T2l[T3l . The motion of carbon- 13 nuclear spins can also be detected 
optically in the NV- system 031, and it has even been possible to map the NV- electron spin 
state onto a nearby C-13 nucleus, and get it back again 031 . In semiconductor quantum 
dots the coherence of electron spins has been optically controlled [fl6l and initialized ifTTTl 
and tunnel coupling between two electron spins in neighbouring quantum dots has been 
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Qubit 


Control 


Method 


Reference 


Electron donor spin 
QD electron spin 
Nuclear donor spin 
Molecular nuclear spin 


Electron donor spin 
Exciton 
Electron spin 
Electron spin 


Optical excitation 
Optical excitation 
Electric Field 
RF pulse 


Ref. 

Ref. a ism 

Ref. a 
Refs. JHHa 



Table 1. Summary of proposals for generating pairwise entanglement of qubits with the help 
of controlled intermediaries. Different approaches will be optimized by choosing different 
material systems. 

detected optically |[T8l . Remote spin coupling through a spin bath has been demonstrated 
in a lithographic quantum dot system Ifl9l . Electrons spin states have been used to manipulate 
nuclear spin qubits in other systems as well [|20ll2Tl . in a time shorter than needed for direct 
addressing. 

In this paper, we consider explictly a system for entangling two such long lived spin 
qubits, via a third, central, electron spin of a different species that might have a much shorter 
decoherence time. This central spin can be optically excited, and in the excited state the 
electron wavefunction typically has a greater spatial extent than in the ground state (see Fig. 
1). If this larger wavefunction overlaps with the two neighbouring spins, it gives rise to an 
exchange coupling. This possibility was first raised by Stoneham et al. Q who introduced 
a scheme for entangling deep donor spins in silicon, where various defects could be used for 
the different spin species (for example. Mg + , Se + , or Bi + |f3l|22l are good candidates). In 
contrast to the previous work, we shall here study a more general situation in which there is 
only limited control over the various coupling parameters in the problem, and show that it is 
still possible to obtain a highly entangling gate operation. 

2. Model 

Let us first introduce a general model for the system we are considering. It has the following 
Hamiltonian, in the usual notation using Pauli spin operators oQ: 

H = E Q af + E c af + E Q ,af + 

|e> (j^aQ ■ a c - aaf ■ a°) + J 2 (o& ■ a c - aaf ■ a c z ) + c^ ) (e|, (1) 

Q and Q' label the two qubit spins; C is the central (control) spin which has two degrees of 
freedom: one is its spin a c and the other is its orbital state which we restrict to the space 
spanned by states {\g), |e)}. We assume an allowed optical transition of energy u between 
\g) and |e), but that Q and Q' do not couple directly to an optical field. Each of the Ej gives 
the Zeeman splitting of spin j in an external magnetic field of strength B {Ej = fijB). J ls 
J 2 is the exchange coupling between spins Q or Q' and C respectively, which is only present 
when C is in the excited state |e). For a = it takes an isotropic Heisenberg form, and for 

| We use a matrices that all have eigenvalues ±1. 
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(a) 




Figure 1. In the optical ground state of the central control qubit (b), the wavefunctions of 
the three species do not overlap and there is no spin-spin exchange interaction. When the 
central control is optically excited (a) its wavefunction has a larger extent and so activates the 
spin-spin coupling. 

a = 1 it represents an XY type coupling. We assume that the control-qubit coupling strength 
when the control is in the \g) state is negligible in comparison with the coupling when the 
control is in the |e) state. Calculations in Ref. Il23l show that this is indeed the case for donors 
in silicon, where the |e) coupling can be two to four orders of magnitude larger. We also 
ignore direct donor-donor coupling, which again is valid for the silicon donor system. For 
suitable spatial configuration with donors separated by around 25 nm, direct donor - donor 
coupling can be as low as 4.73 x 10~ 3 GHz whereas the control mediated coupling strength 
is still up to 157 GHz ||24ll25l 

In order to control the interaction, a laser is applied with frequency oji, and so introduces 
an oscillatory term into the Hamiltonian. The oscillation can be removed by transforming into 
a frame rotating at oji and making a rotating wave approximation, whereupon we write: 

H = E Q af + E c af + E Q ,a?' + ^ (\e)(g\ + \g)(e\) 

+ |e> (ji(a Q ■ a c - aaf ■ a c z ) + J 2 (a Q ' ■ a c - aaf ■ a c z ) + a) (e{2) 

where A = cu — ui is the laser detuning from the transition and is (generally time- 
dependent) Rabi frequency. 

A general state of the three spin system |0) is given by a superposition of the spin 
basis states |{|, |})q <g> |{|, l})c ® |{t> where the arrows represent the spin up or 

down projection along the ^-quantization axis. For convenience, we adopt the usual qubit 
notation \QCQ') with Q, C, Q' being either for the 'down' or 1 for the 'up' projection of the 
respective spin qubit. 

We shall be concerned with the situation in which we initialize the system in |g) and 
where qubit C is prepared in the state |0). We then allow the system to evolve under laser 
excitation until there is a 'revival' such that C returns to the state |0). We will show that the 
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H 3 
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1 1 1 1> 


H 2 


1 


{|110>,|101>,|011>} 


#1 


-1 


{|100>,|010>,|001>} 


H 


-3 


|000) 



Table 2. Table showing the four uncoupled subspace and the notation used for each. 

remaining two qubits Q and Q' can become entangled by such an operation. In the following, 
we will explore two contrasting scenarios, First, a fast 'dynamic' optical excitation in which 
the system is excited suddenly by a pulsed laser, then allowed to evolve for a time before 
sudden de-excitation. Second, we will look at an 'adiabatic' approach in which the laser 
intensity and/or frequency is changed slowly and continuously such that the system follows 
its instantaneous eigenstates. 

3. Dynamic Excitation 

Consider first a laser that is resonant with the \g) — |e) transition (i.e. A = 0). If a rectangular 
pulse is applied for a time t\ — 7r/fi all population is transferred from \g) to |e) thus activating 
the spin couplings. After waiting for a specified amount of time this interaction is deactivated 
again using an identical pulse. For the dynamical approach to work, we require a system 
where a — 1, i.e. that the Ising part is removed from the Heisenberg interaction and we 
are left with an XY coupling. We will discuss the reason for this at the end of this section. 
Hamiltonians not satisfying a = 1 are not amenable to the dynamic method and the more 
general adiabatic approach discuss later must be used. 

Let us assume that the optical excitation is fast in comparison to the subsequent spin 
dynamics, which are described by the restricted Hamiltonian (e\H\e): 

H e = E Q af + Ecaf+E Q ,af+(j 1 (oS ■ a c x + ■ a c y ) + J 2 ( ff ? • a c x + o? ■ a c y )) .(3) 

H e conserves the total spin projection: Y, z = af + erf + af . Therefore the evolution can be 
partitioned into subspaces of different E 2 : 

H e = H © H x © H 2 © H 3 (4) 

where Hi is the Hamiltonian of the subspace with S z = 2i — 3 (see Table 2). The central 
(control) qubit is set to |0) initially and so we need not consider the H 3 space further. Let us 
now analyze the dynamics of the other subspaces. 
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3.1. Hi subspace (H z — —1) 

In the basis of states { 1 010) , 1 100) , |001)}, and for qubits with the same g-factors (Eq = Eq>) 
we have 

( R J[ J' 2 \ 

(5) 



H 1= E 



c 



J 



J[ 1 
\ J 2 o 1 

with R = (2Eq/E c ) — 1, J[ = 2Ji/E c and J' 2 = 2J 2 /E C . There is always one eigenvector 
that is orthogonal to \A) = |010): 

J(|100) - 41001) 



\E) 



(6) 



Let us define a state that is orthogonal to both \A) and \E): 
J' 2 \IW) + J(|001) 



\T) 



J'i 2 + Ji 



1 2 



(7) 



and rewrite the Hamiltonian in the basis {\A), \T), \E)} 

( 

Hi = E c 



R 



l2 



J'" 



J[ 2 + 4 2 



V 



o^ 
o 



(8) 



1 / 





Prior to laser excitation, the system contains no component of \A)\ it can therefore be written 
as a superposition of \E) and \T). \E) is an eigenstate so only accumulates phase. \T) 
undergoes Rabi cycling to \A); after each cycle all population returns to \T), and the system 
'revives' such that no excitation is left on the central qubit. This 'revival time' t rev is given 
by: 

2nix 



tr 



l) 2 + 4J 1 ' 2 + 4jf 



(9) 



E c ^(R-iy + AJ[ 

where n is some integer, the number of oscillations that have occurred. At revival, we have 

\T) ^e lf?T |T), 
\E) ^e i0E \E) 



(10) 



where 



E c {l + R)t T 



Ent 




l) 2 + W 2 + 4J. 



r/2 



C'rev 



(r-i) 2 + aj[ 



l2 



w 2 2 



(11) 
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3.2. H 2 subspace (E 2 = 1) 
In the basis { 1 101), |011), |110)} we can write 

/ -R J[ J' 2 \ 



Ho = Ei 



c 



\ 



4 







-1/ 

After the same revival time t TCV , C again returns to the |0) state, such that the state | A') 
undergoes the following transformation: 



\A>) 



where 



9 a' = nn + 



'\A'), 

E c (l + R)t T 



wn I 1 + 



(R+l) 



(R - l) 2 + AJ[ 2 + AJ> 2 



(12) 
|101) 
(13) 

(14) 



3.3. H subspace (T> z = —3) 
Finally, we have 

1 000) e i6z \000), 

where 

9 Z = -E c (R + 2)t r 



2{R + 2)mr 



{R-lf + 4J{ 2 + AJ' 2 ' 



(15) 
(16) 



3.4. Evolution of the logical qubits; entangling power 

Combining the dynamics for the different subspaces, the overall unitary evolution of the 
system in the logical basis of qubits Q and Q' ({|00), |01), |10), |11)}) is 



U' 



J9 2 









Ai A 2 

A 2 A 3 









e i0 A i ^ 



\ 



(17) 



where 



Ai 
A 2 
A 3 



e idE J[ 2 + e WT J' 2 2 



J[ 2 + 


7/2 

J 2 




- e idE ) 


J[ 2 + 


7/2 

J 2 


e i6T J[ 2 + 


e ie zj> 2 


J[ 2 + 


7/2 

■ J 2 



(18) 



To determine the extent to which this evolution creates entanglement, we use the measure 
of average gate entangling power developed by Zanardi, et al. [|26l . who considered a bipartite 
state l^), which lives in a Hilbert space 7Yi ® 7Y 2 - The entangling power of U is found by 




taking the average of the linear entropy of the reduced density matrix (px = tri[|\&)]) over a 
uniform distribution of input product states ® |^ 2 ): 

e(U) = E{U\^)®\^ 2 )) M \ (19) 

E(\ty)) = 1 — tr(p! 2 ) is the linear entropy of p\. The maximum value of the entangling power 
is about 0.22, and it falls to zero for a gate that produces no entanglement. Using Eq. 5 from 
11261 we determine that 

e(U) = 1(8 - 2|A!| 2 - |A X | 4 - 4|A 2 | 2 - 2|A 2 | 4 

lo 

- 2|A 3 | 2 - |A 3 | 4 - 2Re[e l( ^+^' ) A 2 2 ] 

- 2Re[e i(9z+ ^' ) A 1 A 3 ]). (20) 

Fig. |2]shows the average entangling power of the gate, e(U), after the first revival n — 1, 
for R — 1. [§| The entangling power drops to zero when either Ji = or J 2 = 0: If either 
of the qubits is not interacting with the central control, then no entanglement is possible. By 
contrast e(U) is maximized (and reaches its theoretical maximum) when J[ = J' 2 - and in 
Fig. [3] we show e(U) for different values of the ratio R in this equal coupling case. As R get 
closer to unity, the entangling power is larger for smaller J[ = J' 2 . For larger J[ = J' 2 all plots 
approach maximal entangling power. Overall, we can conclude that the dynamic gate has a 
reasonable entangling power over a wide range of parameter space. 

We can see that it is essential that a = 1 for the dynamic approach to work; if this were 
not the case then the revival times in the Hi and H 2 subspaces would not coincide, invalidating 
the analysis presented here. In order to overcome this restriction we must change our strategy 
and we shall discuss this next. 

§ For even n the average entangling power vanishes, e(U) = 0, whereas for odd n the entanglement revives to 
the same level as for n = 1. 
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4. Adiabatic Excitation 

An alternative method for creating entanglement in our system relies on adiabatic following of 
eigenstates. It can be implemented by slowly modulating the intensity of a laser that is close 
to resonance with the optical transition of the central qubit. Prior to excitation, the system is 
prepared in a superposition of the computational basis states, \QQ') G {|00), 1 01) , 1 10) , |11)}. 
The laser intensity is then varied such that adiabatic following of eigenstates occurs, so if the 
intensity is decreased again population returns to the computational basis. 

With the laser on, each of the eigenstates consists of some superposition of \g) and 
interacting |e) levels, such that the eigenenergies are determined not only by the optical 
coupling but also by the Heisenberg interaction between the spins. Fig. @] shows such an 
eigenspectrum of Hamiltonian (0) as a function of A/Vt. For A/ ft — > oo, the eigenenergies 
tend to the Zeeman split levels comprising the logical basis. The relative spacing between the 
eigenstates changes with laser intensity - i.e. when A/Q approaches zero. Each eigenstate 
therefore acquires a different dynamical phase as a consequence of its time evolution, which 
results in different final phases of the logical states - and thus enables the implementation of 
a controlled phase gate. 

In order to determine which pulse shape and temporal profile is suitable for achieving 
adiabatic following, we shall for the moment neglect the coupling between the spins in the 
excited levels. This gives us eight uncoupled two level systems (2LS), each of which is 
driven independently by the laser. In this case it is straightforward to derive an 'adiabaticity 
condition' which ensures eigenstate following, with suppressed transitions between the 
eigenstates: 

Q(t)A(t)-Q(t)A(t) 
2[A(t) 2 + fi(t) 2 ] 3 / 2 

as in Landau-Zener theory. This condition can be derived by analysing the time-dependent 
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Figure 4. Example of eigenspectrum as a function of the detuning A for typical parameters. 
The eigenstates tending to the computational basis states for 17/ A — > (i. e. to the far right in 
this figure) are colour-coded as follows: |00): green, 1 10) and 1 01 ) : blue, |11): red. 



unitary transformation of a driven 2LS Hamiltonian in the diagonal basis and stipulating that 
the coupling between eigenstates should be small compared to their energetic spacing |[27l . 
For a Gaussian profile of the laser intensity, 

tt(t) = n exp[-(t/r) 2 }, (22) 

and constant detuning A, inequality (l2~TT) can be satisfied by demanding f2o/A 2 <C r. 
Adiabatic following is therefore always achieved in the limit f2 <C A together with a 
sufficiently large pulse duration r. [jj 

The spin-spin interactions mean that the system cannot be regarded as eight separate 2LS. 
Rather, the eigenstates are coupled states which ultimately generate the desired entanglement. 
Nonetheless, the inequality (l2~TT) is still a requirement for achieving adiabatic following, but it 
is not always sufficient; it is also essential that the eigenstates belonging to the computational 
subspace without laser irradiation must be energetically distinguishable from those outside 
this subspace. This avoids population leakage from of the computational basis associated with 
a mixing of eigenstates - and is reasonable in our scheme which presupposes two different 
species for Q and C. 

4.1. Action of the adiabatic operation 

As in Section [3l it suffices to analyze different S z subspaces separately. Unlike for the 
dynamic excitation, however, the restriction a = 1 is not a requirement for the adiabatic 
scheme (and neither is Eq = Eqi). Once more, the control qubit should be initialized to 
a z = —10. The zero excitation subspace then only contains the logical 1 00) state. Similarly, 

|| For small 17o/A, t must be automatically large since the interacting levels are only weakly excited. 

It is important to have the central qubit in a well-defined initial state because of subspace dependent energy 
shifts of the eigenstates. In general, this leads to different pulse durations for a successful entangling operation 
depending on whether C is initially in the a z = — 1 or the a z = +1 state. 
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the two excitation subspace contains only the logical |11) state. Therefore, no population 
transfer between these and other logical states is possible and each merely accumulates a 
phase after the adiabatic pulse has been applied. The situation is different for the single 
excitation subspace, which is populated by the two interacting states 1 01) and |10). For 
this subspace, a more complex unitary operation between the logical states results from the 
adiabatic operation. 

The most general action of the adiabatic operation can be therefore described by the 
following unitary matrix in the basis of the logical states {|00, |01), |10), |11))}: 



/ g*</>00 



u, 



ad 



\ 











i> 

X 






x' 









3*011 



\ 



/ 



(23) 



where 4>ij is the phase acquired by the state The coefficients ij), ip', X an d x' form a 
unitary 2x2 submatrix accounting for population transfer between 1 10) and 1 01) as well as 
the phase acquired by each of these two states. The complexity of the Hamiltonian © makes 
it difficult to find analytical expressions for the elements of Eq. (1231) . Nevertheless, for a given 
set of system and laser control parameters U a( i is straightforward to obtain numerically. 

The structure of U a d takes the same form as Eq. (flTI) obtained in Sec. [3] for the dynamic 
operation. This enables a direct comparison of the entangling power of the dynamic and 
adiabatic approach using the measure defined in Eq. (1201) . 



4.2. Adiabatic CPHASE Gate 

Under certain conditions the off-diagonal terms in Eq. ([23]) are zero. Using the numerical 
techniques discussed earlier, we find that this is the case where either: 

(i) there are degenerate logical qubits Eq = Eq/ with equal coupling J\ = J 2 to the control 
qubit, or 

(ii) there are non-degenerate logical qubits Eq ^ Eq/ and a pulse duration r longer than 

(Eq - Eq,)- 1 . 

In this case the adiabatic operation is simply 



/ 


g«0OO 











\ 







g«0Ol 


















gi<AiO 







V 











g«011 


/ 



(24) 



which is locally equivalent to the CPHASE gate when B1I281: 

V? = 000 - 001 - 010 + 011 = 7T. (25) 

This condition that can always be satisfied by choosing an appropriate pulse duration. 

We now explain why a system which satisfies the less restrictive set of conditions (ii) 
above gives a unitary operation of CPHASE form. Our explanation need only consider the 
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i = 1 subspace since the others are in CPHASE form under any conditions. We write a 
general initial state in this subspace characterized by the amplitudes a and (5 as 

|V(0)) = a|100)+/3|001). (26) 

Adiabatic following of eigenstates means that this state evolves under the influence of the 
laser into 

m)) = aMt))+PW(t)), (27) 

where \fJ>(t)) tends to 1 100) and \u(t)) to |001) as the Rabi frequency goes to zero. The 
time evolution of the slowly changing constituent eigenstates of Eq. (1271) follows 

|/i(t)) = e^ i |/i), (28) 
\u(t)) = e iE » t \v), (29) 

where fx and v denote the instantaneous eigenstates and and E v are their associated 
eigenvalues. The time evolution of \ip(t)) may thus be written as 

= e lE ^ (a\fi) + e i{E »- E ^[3\v)) . (30) 

The i = 1 subspace consists of three spin states in each of the ground and excited optical 
states, so that \p) and \v) will be composed of up to six different states. Focussing on the 
physical states which correspond to the two logical states, \p) and \v) can be written as 
follows: 

|^> = (p|100> + #01))<g)|<7> + |m>, (31) 
\v) = (r|100) + s|001» <g> \g) + \n), (32) 

where p,r,q and s are appropriate amplitudes of this decomposition and \m) and \n) 
contain the contributions of the four remaining states |010)|g), 1 100) |e) , |001)|e)and 1 010) | e) . 
Inserting Eqs. (|3~Tll32l into Eq. © yields 

= e lE ^ (a\m) + e i{E »- E ^ (3\n) 
+ (ap + e^-^pr) |100) 

+ (aq + e i{E "- E ^ps) |001) ) . (33) 

At the end of the pulse, adiabaticity ensures that the first two terms of Eq. (1331) disappear. The 
third and fourth terms show cycles of constructive and destructive interference of the states 
1 100) | g) and |001)|g). The interference oscillations range between ap ± (3r for |100)|g) and 
aq ± (5s for |001)|g) with a period of [E v — E fJ )/2iv = AE/2n. Our simulations show that if 
the interference period AE/2% is considerably faster than the pulse duration r, all population 
is restored to the original levels at the end of the pulse, as illustrated in the scenario of Fig. 
|5] - which of course allows the realisation of a CPHASE gate with suitable laser control 
parameters. Conversely, comparatively fast laser pulses generally transfer population between 
the two logical states, leading once more to the gate described by Eq. (|23T) . 
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Figure 5. Interference of population between |10)l and 01} i in the i = 1 subspace. This 
effect is a consequence of the time evolution of a superposition of eigenstates as explained 
in the main text. As shown, a sufficiently long pulse duration damps the oscillations out and 
restores all population back into the original levels at the end of the pulse. The Gaussian pulse 
is centred around t = 0, where time is given in units of the pulse width r, which in this case is 
set to 150 ps. 



4.3. Entangling Power 

We simulate the adiabatic operation by integrating Hamiltonian Eq. <0Q) with a Gaussian 
profile of the Rabi freguency as in Eq. (|22j) . In order to prevent Landau-Zener transitions 
between the eigenstates, r needs to be suitably large. Depending on whether a pure phase 
gate or a more general entangling gate is desired, the effect of the operation can be obtained 
after the pulse has finished by extracting either a non-trivial phase as in Eq. (l25l) or the 
unitary matrix Eq. (T23T) . The average entangling power of both these quantities can then be 
determined using Eq. (|2Q|) . 

We find that a more pronounced difference between the Zeeman splittings of all three 
spins makes the adiabatic following more robust and permits the application of a pulse with 
shorter duration. This might be achieved by using species with varying g values, or through 
an inhomogeneous magnetic field. 

Fig. [6] shows a typical plot of the average adiabatic entangling power (as in Eq. 
( fl9l )) as a function of Ji and J 2 , and Fig. [7] presents a cross section along the along the 
diagonal J\ = J 2 . El As for the dynamic gate, the most entangling region occurs along 
the diagonal where J\ and J 2 are equal; in contrast to the dynamic gate, the graph shows a 
rather complicated oscillatory structure. This is connected to how phase accumulates in the 
adiabatic gate - which is in turn related to the values of J\ and J 2 . The plots contain some 
areas where adiabatic following does not occur, causing significant population leakage away 
from the computational basis states. In this case e(U) is ill defined. However, a well defined 
adiabatic operation can always be recovered by making adjustments to the choice of laser 

+ In order to compare these simulations to those for the dynamic gate we choose more restrictive conditions than 
are necessary: Matching onsite energies Eq = Eq> =0.1 ps -1 , XY type coupling (a = 1) and Eq = 0.1 ps -1 . 
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R=1.2 




Figure 6. Average entangling power e(U) of the adiabatic two qubit gate. The white areas of 
the contour plot (light blue patches without grid lines in the 3D plot) correspond to parameter 
combinations for which the adiabatic gate leads to population leakage out of the computational 
basis, making e(U) ill-defined. J\ and J2 are given in units of 0.1 ps _1 . The pulse duration 
r = 0.5 ns, the detuning A = 0.5 ps" 1 and coupling strength fig = 0.3 ps _1 . N. B. For a 
particular Ji, J2, it is possible to optimize the speed of the gate by varying r, A and ft 




Figure 7. Average entangling power of the adiabatic gate when J\ = J2 (in units of 0.1 ps x ), 
for different values of R. 



control parameters. 
5. Decoherence 

We shall now discuss the effect of decoherence on our predictions. The dominant decoherence 
source or sources will be different for each physical implementation of our Hamiltonian Eq.[D 
and a full discussion of each possible process is beyond the scope of a single paper. However, 
we shall discuss decoherence that arises from spontaneous decay of the optically excited state. 
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Figure 8. Comparison of final population in the computational basis (squares) and final purity 
of the system's density matrix (circles) for the adiabatic gate (orange) and the dynamic gate 
(blue). The decay rate is given in units of ns _1 and the system parameters are J\ = J2 — 
0.05 ps _1 and R = 1.2. Typical laser control parameters have been used, such that both gates 
achieve an entangling power very close to the maximal value of e(U) ~ 0.22. 



This could either be radiative or non radiative If29ll30l . and would be the dominant source for 
deep donors in silicon Il3ll22l. 

We use a standard quantum optical master equation OT1 to model the decay affecting the 
control qubit Q 

P = ~i[H, p] + T (o-po + - ^(a + a^p + pa+a^ , (34) 

where p is the system's density matrix, r is the decay rate (the inverse of the natural lifetime) 
and o + and ct_ are raising and lowering operators with respect to \e) and \g). 

We will explore two figures of merit: the amount of population returned to the desired 
computational basis states, and the purity of the density matrix, after application of the gate. 
For both gate types we have performed a full numerical simulation that for the dynamic gate 
includes the two (rectangular) laser tt pulses. Fig. [8] shows that the fast dynamic gate retains 
a higher purity even for fast decay rates. However, the adiabatic gate is more robust to loss of 
population from the computational basis. 

In Fig. |9l we analyse the dependence of the two figures of merit on J\ = J 2 . For the 
dynamic gate, we keep Q = 5 ps~ x ~ 3 meV constant; this introduces an increasing intrinsic 
error as t TCV decreases and the transient regimes of excitation and de-excitation become more 
important. On the other hand, the purity improves as t rcv becomes shorter when radiative 
decay is included, as shown by the red curves. For the adiabatic gate, a pulse duration 
r = 0.5 ns and coupling strength f2 = 0.066 meV are used, with the detuning adjusted in 
the range A = 0.16 — 0.6 meV to give a maximum entangling power for each of the J\ = J 2 
points shown. As can be seen in Fig. [9l the adiabaticity condition is very well satisfied and 
there is no population leakage in absence of radiative decay. However, decay events inevitably 
lead to some population leakage when included in the model. Unsurprisingly, a stronger XY 
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Figure 9. Intrinsic gate errors (blue curves) and additional effects of decoherence (red curves): 
the final purity and the final population of the computional basis states are shown as a function 
of the coupling strength J\ = J2 (in units of 0.1 ps _1 ). A decay rate of 1 ns _1 , R = 1.2 and 
typical laser pulse parameters have been used. 



interaction improves performance for both purity and loss of population; this is in contrast 
to the dynamic scheme, for which the loss of population gets worse when the interaction 
strengths are larger. 

6. Conclusion 

We have shown that it is possible for a central control to be a mediator of entanglement 
between two qubits, even if the coupling strengths between mediator and qubits are different. 
When a dynamic approach is taken, the coupling must be of XY form - but if an alternative 
adiabatic gate is performed any coupling form is permissible. Further, the gate is close to 
maximally entangling over a wide range of parameter space, in both cases. 

The proposed protocol could immediately be tested in a range of experimental systems, 
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include molecules with coupled electron and nuclear spins, and donors in silicon. Possible 
experiments would include the demonstration of remote entanglement between two centres 
that are not directly coupled, or the demonstration of a simple algorithm. 
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